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LATTICE POINTS ON CIRCLES AND DISCRETE VELOCITY 
g ; MODELS FOR THE BOLTZMANN EQUATION 

o : 

<n . laura fainsilber, par kurlberg, and bernt wennberg 

1. Introduction 

The phase space density / of a dilute gas evolves according to the Boltz- 
mann equation. In the physically relevant case, the gas would be confined to 
■ a subset f2 C IR 3 , and then f(x, v , t) : x IR 3 x IR+ — > R+, where x denotes 

a position in space, v G IR 3 is a velocity, and t denotes the time. From a 
^ \ mathematical point of view, it is equally natural to consider the Boltzmann 

equation in any spatial dimension, and in some cases because of symmetries 
of fi, it is also relevant to consider f2 C IR dl and v G IR d2 with d\ < d 2 . 

By a dilute gas we mean one where the particles interact with each other 
essentially only by pairwise interactions. Moreover, the Boltzmann equa- 
tion assumes that the particles are so small compared to other distances, that 
O- ■ they can be considered to be points. 

! Under these hypothesis, one can formally derive the Boltzmann equa- 

O ; tion (see 0) 

(1) d t f(x,v,t)+vVJ(x,v,t) = Q(f,f)(x,v,t). 

^ ■ The left hand side describes the evolution of the density by free transport, 

and the right hand side describes the impact of collisions. Per definition, a 
collision is a pairwise interaction that takes place instantaneously and at one 
single point in space. Hence x and t appear only as parameters in Q(f, f), 
and we can write 

Q(f,f)(v) =11 (f(v')f(vl)-f(v)f(v*))q(\w\,co S e)dS(u)dv* 
(2) 

where the velocities "before and after a collision" are related by 

v' = \{v + v*) + \w\ u 

(3) 

K = \{v + f*) — \w\u, 
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with w = (i>* — v)/2, and with cos6> = dv* is the Lebesgue measure 

in R. d , and dS(u) is the surface measure on S^ -1 . Note that the pair of 
velocities before a collision, v and t>*, and the pair of velocities after the 
collision, v' and i>£, are the endpoints of a diameter on the sphere which has 
its centre at and diameter |u* — v\ . This is exactly the condition needed 
in order that the collisions preserve the momentum and energy of the pair 
of particles. For d — 2, the sphere becomes a circle, and this motivates the 
title of the paper. 

In a discrete velocity model (DVM), the velocities are concentrated on a 
(usually finite) set of points Vj G M d in the velocity space: 

f(x,v,t) = ^2fj(x,t)6 v=Vj . 
j 

The Boltzmann equation ^ is then changed into a nonlinear system of 
conservation laws, 

(4) d t fj + Vj ■ VJ, = T ik' iUh' ~ fih) > 

k,k',f 

where the constants rj £ > must be chosen so that (0|) makes sense from 
a physical point of view. In particular we require that (vj, Vk) and (vy, v^) 
define two diameters on the same sphere, just as for the usual Boltzmann 
equation. 

The first example of a discrete velocity model is that of Carleman (|4|), 
which has two velocities in K. Many other models have been proposed, 
and there is a large literature on how to construct and analyse physically 
realistic models (i.e., that satisfy the right conservation laws and an entropy 
principle), see eg.0. 

Besides offering many interesting mathematical challenges (for exam- 
ple, there is no general theory of global existence of solutions to systems 
like @) the DVM:s are also candidates for the numerical approximation 
of the real Boltzmann equation O. This leads naturally to the following 
question, which is the subject matter of the paper: 

Suppose that we choose the discrete set of velocities to be hZ d , i.e. the 
integer lattice in IR d , scaled by a factor h, and that we take 

f h (v) = h,hS v =H > 

so that f h — > /, in some suitable sense, where / G Z/^R^). Is it then true 
that Q(f h , f h )(v) -> Q(f, f)(v) for all v G hZ d when h -> 0? 

This property, which is called consistency, together with stability are 
main ingredients when proving that a numerical method converges. 
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The answer is yes. This was proven by Bobylev, Palczewski and Schnei- 
der for dimensions d > 3. In this paper, we prove that it is also true 
for d = 2, and hence for all relevant cases. 

Results of this kind are interesting, because they provide examples that 
are relevant to previous results of Desvillettes and Mischler ([ 9 1), who proved 
that solutions to families of DVM:s can converge to DiPerna-Lions' solu- 
tions to © if certain conditions are satisfied. 

Our result should not, however, be considered as relevant for numerical 
analysis, because the rate of convergence is so slow that a numerical method 
based on the theory presented here would hardly ever become useful. 

The family of models considered here can be seen as coming from a 
rather straightforward discretization of the collision integral ©. This inte- 
gral should be interpreted as an average over the 2d — 1 -dimensional mani- 
fold defined by 

(5) M v = {(v*,v',v' t ) e R 3d s.t. v' + vi-v* = v 

I /|2 i I / |2 I |2 I 1 2 1 

and © is an iterated integral over this manifold. For a fixed v, we write 
w = (t>* — v)/2, and then ® becomes 

v' = v + w + \w\u 
v'x = v + w — \w\u 



and also v* = v + 2w. We then write 

(6) g v (w,u) = (f(v')f(v'J-f(v)f(v*))q(\w\,cos6), 

and so (after changing variables in the integral), 

Q(fJ)(v) = 2 d [ ([ g v (w,u)dS(u)) dw. 

If g is sufficiently regular (continuous), and decays sufficiently rapidly for 
large w, then the Riemann sum for the outer integral converges: 



(2h) d T I 9v(hC,u)dS(u) 



(7) — > 2 d / ( / g v (w,u)dS(u) ) dw 

Jr n KJs' 1 - 1 ' J 

when h — > 0. In order to construct a consistent DVM, it is then sufficient 
to evaluate the inner integral in terms of the values of g on the lattice points 
hZ d , in such a way that the result converges to J sd -i g(w, u) dS(u). While 
with the formula ©, the collision integral should be taken over all u G 
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S d \ we have here only access to those u for which v' and belong to 
KL d . But this is automatically achieved if £ £ Z d , and if u = C'/\C% where 
C £ Z d and |C'| = |C|; then for all v £ hZ d , 

v + h(±h\(\u £ hZ d . 

However, note that with this construction, the center of the sphere is re- 
stricted to lie on a lattice point, and so it excludes cases like v = (0, 0), 

v* = (h, h). 

Giving all points on the sphere equal weight, one arrives at the expression 

( g ) E (/(«W.) - /(«)/(«•)) q(\K\,cose), 

Td[l ^ } c& d 
IC'NICI 

for approximating the inner integral in (JT]). The function rd(n) denotes the 
number of points with integer coordinates on a sphere in W 1 with center at 
the origin and radius y/n, i.e. the number of integer solutions to x\ + ■ • ■ + 

x 2 = n. 

We write, for all v £ hZ d . 

Q h (f,f)(v) = 

< 9 > (^Et?^ E (f(vVM)-mf(v*))q(\h(lcose). 

IC'NICI 

In the two-dimensional case, all the terms in the sum are 27r-periodic 
functions of 6, and assuming sufficient regularity, they can be expressed as 
a convergent Fourier series. It is then natural to introduce the exponential 
sum 

(10) S(n,k)= eik<>u 

udlP-:\u\ 2 =n 

where 9 U is defined by u = \u\ ■ (siad u , cos6 u ). We will see in Section|4| 
that to prove that dU) converges to the angular integral in Q, it is enough to 
prove that for k ^ 0, the terms S(n, k) converge to zero sufficiently fast as 
n — ► oo. Similar exponential sums are relevant for any dimension, and the 
work of Bobylev et al. also involves such estimates. 

Here the needed estimate is given as Proposition|^lin Section|3J Then in 
Section|Uwe put the estimates togheter to a proof of the main result: 

Theorem 1. Consider the Boltzmann equation in two dimensions. Assume 
that f and q are so smooth that the function g v (w,u) defined in © is a 
C 2 -function. Then for all v £ KL 2 

\Q(fJ)(v)-Q h (fJ)(v)\^o 
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when h — > 0. 

Section finally, contains a numerical illustration on the distribution 
terms where the circles contain many points, and we indicate how the com- 
putational cost could be reduced without loosing accuracy (still without the 
intention of actually giving an effective algorithm). 

A more general construction of discrete velocity models on scaled inte- 
ger lattices hi 2 consists in finding sets of integer points on the manifold 
M. defined in ©. In this way, mass and energy conservation are automati- 
cally satisfied, but one also needs to verify that these are the only conserved 
quantities. And finally, in order that the models converge to the continuous 
model when h — > 0, it is necessary that the integer points are more or less 
uniformly distributed on M. . 

The models studied here are constructed by discretizing, one at a time, 
the iterated integrals ©. An alternative way of writing this integral was 
introduced by Carleman \4\. Using that v' — v and v [ — v are orthogonal, 
one can write (here we specialize to d = 3) 

Q(fJ)(v) = 

! I (f(v')f(vi)-f(v)f(v # ))q(w,co S 9)—!—idE(vi)dv', 
Jk 3 Je ViV , \v-v'\ 2 

where E v y is the plane that contains v and is orthogonal to v'—v, and where 
dE(v'^) is the Euclidean measure on this plane. Heintz and Panferov lfl~8l 
have analysed a DVM based on this interated integral, and proved that the 
method is consistent with the continuous model. This is somehow easier, 
because on all planes, the integer points are uniformly distributed, and they 
are all found by solving linear Diophantine equations. However, the density 
of points depends strongly on v' — v, and so it is far from trival to prove the 
consistency. And again, the two-dimensional situation is more difficult, and 
has not yet been studied. 

Yet another approach was introduced by Rogier and Schneider who 
used the theory of Farey series to discretize the angular variable in the col- 
lision integral. 

2. Number theoretic background 

2.1. Points on spheres; Asymptotics. To prove that © converges to the 
correct limit when h — > 0, one has to study the set 

{C/ICI : C G Z d , ICI 2 = "} 

and to show that the points of this set are sufficiently well distributed on 
S^ 1 when n is large; it is here that the number theoretical issues enter the 
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game. Indeed, we can view the set of points with integer coordinates on a 
sphere of squared radius n centered at the origin, 

d 

{(x 1 ,...x d ) e z d ,J2 x2 i = n ^ 
i=i 

as the solution set for a quadratic form, and use the theory of integral qua- 
dratic forms to get estimates on the number of points (see for instance £\3i ). 
The expected number of points with integer coordinates on a sphere clearly 
depends on the dimension d. The naive approach to find the order of mag- 
nitude for a given dimension is to use the volume of a ball, divided by 
the number of spheres contained in the ball. The volume of a ball of ra- 
dius y/n grows as n d l 2 whilst the number of spheres is n. For d — 2, this 
leads us to expect a constant number of lattice points on circles, for d = 3 
a growth proportional to \fn, etc. However, for small d this approach is 
misleading; the growth is quite irregular, and depends on the divisor struc- 
ture of n. For d = 2, we will see below that only values of n of the form 
n = 2 s q 2 p" 1 . . . p" r , where q is a product of primes of the form Ak + 3 and 
the p^s are primes of the form Ak + 1 (see below), yield circles with lattice 
points, and thus most circles have no points at all. In fact, Landau proved 
in 1908 that the number of circles with at least one lattice point, of integer 
squared radius smaller than x, grows as Cx/y/logx. Moreover, there are 
also infinite families of circles with very few lattice points; radii that are a 
power of 2 yield 4 points for instance, and radii that are the square root of 
a prime of the form p = Ak + 1 yield exactly 8 points. On the other hand, 
the number of lattice points on a circle is not bounded, for instance a circle 
with n = pi . . . p r as above where all the pi are distinct from each other has 
4 • T points. 

In dimension 3, all values of n not of the form n = A s (8k + 7) yield 
spheres containing points with integer coordinates. This still leave a fairly 
large number of spheres with no points on, but for our purposes this does 
not really matter, as such spheres do not appear in the summation formulas 
(there is no relevant value for £.) Among the spheres with lattice points, 
multiplying the radius by a power of 2 does not increase the number of 
points, but if we correct for this fact, the ratio between the number of points 
and the naive estimate is bounded, up to constants only depending on e, 
from above by n £ , and below by n~ e for all e > (see |[T3l Ch. 4 for exact 
formulas involving class numbers or L-series.) 

The higher-dimensional cases behave in a somewhat more regular fash- 
ion. Lagrange proved that every positive integer can be written as the sum 
of four squares, and thus for dimension d > 4, every sphere whose squared 
radius is an integer, has lattice points. For d = A the number of points still 
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oscillates rather wildly, with spheres with radius a power of 2 having just 
24 points, but for greater dimensions, the naive estimate gives the correct 
asymptotic growth of the number of points. 

Getting circles (or spheres) with "sufficiently many" lattice points, how- 
ever, is not quite enough for our purposes: we also need that the lattice 
points are sufficiently uniformly distributed when projected on the unit 
sphere. In dimensions 3 and higher, this follows from estimates on Fourier 
coefficients of modular forms. The case d > 4, with some restrictions 
on the set of numbers in which n tends to infinity when d = 4, is due 
to Pommerenke [33. For d = 3, Duke [11| and Golubeva-Fomenko iflU 
used Iwaniec's (TT^l estimates on Fourier coefficients of half integral weight 
forms to obtain uniform distribution. Unfortunately, these techniques do 
not apply in dimension 2. Moreover, there are circles with large number of 
lattice points that are poorly distributed: 

Theorem 2. (Cilleruelo |6|) For any e > and for any integer k, there 
exists a circle x 2 + y 2 = n with more than k lattice points such that all the 
lattice points are on the arcs y/ne^ 2 ^ t+ ^ 1 with \9\ < e, t E {0, 1, 2, 3}. 

On the other hand, we may use some other techniques from analytic num- 
ber theory to show that lattice points on circles are equidistributed on aver- 
age, and this is good enough for our purpose. 

2.2. From points on circles to Gaussian integers. In the plane, we can 
view lattice points on a circle of radius ^/n, centered at the origin, as com- 
plex numbers with integer real and imaginary parts, and squared modulus n. 
It might seem as a trivial restatement, but doing so allows us to use use some 
techniques from algebraic number theory. The Gaussian integers, i.e., the 
set 

Z[i] = {x + iy e C,{x,y) e Z 2 }, 

is the ring of integers of the field Q(i). It shares an important property 
with the ordinary integers, namely unique factorization 1 , i.e., just as every 
integer in Z factors into prime numbers, and the factorization is unique up 
to ordering the primes and multiplying by —1, Gaussian integers factor into 
Gaussian primes, uniquely up to ordering and multiplication by —1, i, —i 
(these and 1 are the units, i.e. the elements having a multiplicative inverse 
in Z[i]). For a more thorough introduction to primes in quadratic number 
fields, see for instance ifTTl . Ch. XV. 

The Gaussian primes (i.e. the elements of Z[i] that cannot be written as a 
product of Gaussian integers with smaller modulus), are of three types: 

'This is rather unusual, the ring of integers in most number fields will not have this 
property. 
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• the prime numbers q G Z such that q = 3 mod 4 remain prime in 
Z[i](e.g. 3,7, 11,19,...); 

• for prime numbers p E Z such that p = 1 mod 4, there exists 
x, y E Z s.t. p = x 2 + y 2 . Hence p factors in Z[i] as a product of 
two Gaussian primes 

p = (x + iy)(x — iy) 

(e.g. 5 factors into (2 + i)(2 - i) in Z[i]) 

• last (and least!), 1 + i is prime (note that (1 + i)(l — i) =2 and that 
1 — i — — i(l + i) is merely "another form of the same prime" just 
as 3 and —3 represent the same prime). 

If n is the sum of two squares, then it can be factored in Z[i] : 

n = X 2 + Y 2 = (X + iY){X - iY). 

If z = x + iy is a prime factor of X + iY, then z = x — iy must be a 
prime factor of X — iY. It follows that prime factors q = 3 mod 4 of n 
must appear in even powers. In addition, multiplying n by an even power 
of a prime g that is congruent with 3 mod 4 changes neither the number of 
solutions to n = X 2 + Y 2 nor the distribution of arguments of the solutions. 

Suppose now that n contains a factor p a , where p = 1 mod 4. The num- 
ber p can be factored in Z[i] as (x + iy) (x — iy), and hence the multiplicity 
of x + iy as a factor of n is a, and the same is true for x — iy. It follows that 
the multiplicity of x + iy in X + iY can be any integer j, with < j < a, 
and the multiplicty of x — iy is then a — j. 

The same calculation can be done for powers of 2; however, the solutions 
given by different choices of j in that case differ by a multiplication by a 
power of i, and so the power of 2 does not influence the number of solutions. 

All solutions to n = X 2 + Y 2 can now be expressed as X + iY = 
y/nexp(i9), where all possible values of the argument 9 can be computed 
as sums of terms deriving from the different factors of n in the following 
way: 

(1) X + iY can be multiplied by any unit, i.e. by ±1 or ±i. This gives 
a term kn/2 in the argument, k — 0, 1, 2, 3. 

(2) If the multiplicity of 2 in n is odd, then the argument must contain 
7r/4, the argument of 1+i; the number of solutions does not change. 

(3) For each prime factor p = 1 mod 4 in n, let a p be the multiplicity 
of p in n, let p = x 2 + y^, and set # p = arg(x p + iy p ). For a 
particular choice of j, < j < a p , the argument added to X + iY 
is jflp - (a p - j)6» p = (2j - a p )9 p . 
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Since the choices of k, and of the different j's are independent, the number 
of different solutions is 4 (a p + l). 



ilmod4 



2.3. Results on the distribution of primes and on the angular distribu- 
tion of points. We will need the following results: 

Theorem 3 (Merten's Theorem, see [ 17 1, Ch. 22.8). 



'[ [ (l-l/p)~e-71og 



x . 



p<x 
p prime 

where 7 ~ 0.57/5 Euler's constant. 

As for the angular distribution of Gaussian primes, a result by Kubilyus 
gives that the angles {0 p } p =i mod 4 are equidistributed in [0, tc/A] in the fol- 
lowing sense: 

Theorem 4 (Kubilyus, ll2"T10 . The number of Gaussian primes uj in the sec- 
tor < a < &rg(u) < (3 < 2tt, \uj\ 2 < u is equal to 



— {(3 — a) [ — — + O (uexp(-b\/\ogi 
7r J 2 logv V 

where b is an absolute positive constant. 

From Kubilyus' Theorem, it is straightforward to deduce (see lfT2ll . p. 92): 



Corollary 5. Ifk e 4N and log A; < by/log x, then 

Y |CQS(A; ^ )I < -logloga; + (l-2/7r)loglogA; + 0(l). 

p=l mod 4 



3. EQUIDISTRIBUTION OF LATTICE POINTS ON CIRCLES 

What is needed for the proof of consistence of the discrete velocity model 
are estimates on the equidistribution of lattice points on circles. 

The aim of this section is to show that lattice points on circles are equidis- 
tributed on average in the sense that the exponential sums S(m, k) converge 
to zero when m goes to infinity. We recall the definition of S(m, k): 

S( m , k) = ^ ■ 

\w'\ 2 =m 
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Proposition 6. If 4 \ k then \S(m, k)\ =0. If 4\k and k ^ 0, there exist C 
and b > such that 



log(^ J2 \ S ( m > k )\) < C - (1 - 2/tt) log 

\ m<X / 

for X sufficiently large and log \k\ < b^logX. 



(log |A:|) 2 



Remark: The mean discrepancy of the angles of Gaussian integers were 
studied by Katai and Kornyei in [20], and by Erdos and Hall in lfl2l . Our 
method is similar to theirs, except that they bound 

1 \- \S(m,k)\ 
X/VhgX^ x r{m) 

instead of 

i^|5(m,fc)|. 

The proof is based on the observation that | S(m, k) |/4 is a multiplicative 
function, i.e. a function / : N — ► C such that f(mn) = f(m)f(n) for all 
m, n such that (m, n) = 1. It turns out that the mean value of a multiplica- 
tive function, under fairly general circumstances, can be bounded in terms 
of an exponential of a sum over primes. To make the paper more self con- 
tained, we include a weak form of the Halberstam-Richert inequality (cf. 

CI). 

Theorem 7. Let f be a nonnegative multiplicative function such that 
(11) £/(tO = 0(o;), 

n<x 

and f{p k ) = 0{k) for all primes p and k > 1. Then there exists C > 
such that 

I^ / ( ro) < C .exp(^i^) + 0( 

m<X \P<X ^ / 

for all sufficiently large X. 
Proof Following Wirsing ll231l . let 



logX' 



n<t 



Then 



[ X 1 ~Y L dt = F(X) logX + 0(1) - f( n ) l °Sn. 

n<X 
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On the other hand, by assumption, we have F(t) = 0(t), thus 



j* E^± dt = (x) 



and hence 

F(X) \ogX < 0(1) + X + /( n ) lo § n • 



n<X 

2 . 



Using logn = Yld\n A(rf)> wnere A is the von Mangoldt function we have 



^/(n)togn=X)/(n)E A ( d ) = E A ( d ) E ^( dm ) 

n<X n<X d|n d<X m<X/d 



m) 



(12) = 5>(cQ ^ /(dm) + J>(d) £ /(<* 

(i<X m<X/d, d<X m<X/d, 

(m,d)=l (m,d)>l 

Now, since A(c?) = unless d is a prime power, we have 

(13) £)A(d) ^ /(dm)= ]T log(p) ]T /(p fe+ 'm) 

ci<X m<X/d, p k+l <X m<X/p k+l 

(m,d)>l k,l>l (p,m)=l 

= i°s(p)f(p k+l ) E /( 

p k+l <X m<X/p k+l 
k,l>l (p,m)=l 

By the assumptions on /, 



/(p**) £ /(m)<O(* + £ /( m ) = 0^ + i)^V 

m<X/p k+l m<X/p k+l ' 

(p,m)=l 

and thus the second term in (fT2"l) is 

= 0(.Y) , 



o \ ^sipW— 

p n <X p 
n>2 



since 

log(p) 



£ E log(p)n 2 p- B < E ^ E ( 2 + -) 22 " m < 00 • 



p n>2 p m>0 



That is, A(d) = logp if d = p fe and fc > 1, otherwise A(d) = 0. 



mi 
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As for the first term in (fl2h . we have (recall that / is multiplicative and 
nonnegative) 

£A(d) f(dm) = J2Md)f(d) E /( 

d<X m<X/d, d<X m<X/d, 

(m,d)=l (m,d)=l 

< E /h E A ( rf )/( rf ) • 

m<X d<X/m 

Now, 

E A ( rf )/( rf ) = E Mp)/(p*) < E MpW) = (*/ 

d<X/m p k <X/m p k <X/m 

k>l k>\ 

since 

E log(p) = 0(X/m) 

p<X/m 

by the Prime number theorem, and 

E k log(p) = O ((X/m) 1/2 log 3 (X/m)) = 0(X/m) 

Thus, 



p k <X/m 
k>2 



J2fn E A(rf)/w = o(E/( m )^) • 

m<X d<X/m \m<X J 

But since / is nonnegative and multiplicative, we have 

E^r^ U(i + f(p)/p + f(p 2 )/p 2 + ---) 

m<X p<X 

< n (c 1 + /w/p) ■ ( x + /^w + /(pw +•••)) 

p<X 

and since 

E (/(pW + /(pW + •••)< E E ^ < - • 

p<X p fc>2 ^ 

we find that 

E?=°(nc+/w/f))- 

Thus, 



F{X) logX = O ( X + X ■ H (1 + f(p)/p) j 
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hence 



X 



o 



i , n P < x (i+f(p)/p)\ 



logX 



logX 



Now, by Merten's theorem, we have 

ri( i - i /p) 



and thus 



F(X) 
X 



o 



p<X 



logX 



n i 



p<X 



O 



(i^ + exP fe 



log A" 

| /(P)~l Hp) 
P 

f(p) ~ 1 



V 



□ 



Proof of Proposition® To see that \S(m, k) /4| is a multiplicative function, 
it is enough to recall the factorization of m into Gaussian primes. Namely, 
if p" 1 , ...p°j J are all prime factors of m with p = 1 mod 4, 

3 qi aj 

S ^ ^ _ . . . e ife(eo+K-2ii)0 P1 +...+( aj -2j. 7 )e PJ ) _ 

<?=o ji=i ij=i 

Here #o is a multiple of 7r/4 which comes from powers of 2 in m, and 
the 9 p . can be computed from the Gaussian factorization as described in 



Section |2~21 Also, because ^ 



\S(m,k)\ 



3 M 



4 if 4 A; and zero otherwise, 



. . . e ik((ai-2h)e n +...+(aj-2jj)8 Pj ) 

and this sum clearly factorizes, each factor containing a sum of terms cor- 
responding to one of the prime factors p. Hence 

is a nonnegative multiplicative function, as stated. In addition it satisfies 
fk( m ) < r(m)/4 for all m. Thus, since 

^~V(n) = e Z : x 2 + y 2 < T}\ ~ vr (v 7 ^ = vrT 



we have 



OCT). 
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Moreover, if p = 3 mod 4 then 



(14) Hp 1 ) 

and if p = 1 mod 4 then 



1 if I is even, 
if I is odd, 



(15) 



hip 1 



ik(l-2j)6 p 



and thus /fc(p') < / + 1 for all prime p and I > 1. The assumptions in in 
Theorem|7]are thus satisfied, and we obtain 



l - \S(m,k)\ = ± £ /fc(m ) < Cexp (j2 

m<X m<X \P<X 



X 

Now, by (fl4l and dl3T >. we have 

/*(p) 

Hence 



fk(j>)-l\ +0 ( I 



P 



loeX 



2 1 cos(k9 p ) | if p = 1 mod 4, 
if p = 3 mod 4 . 



E 

By Corollary 



fkjp) ~ 1 
P 



E 



2\cos(k9 p )\ 



p<X 
p=l mod 4 



P 



E 1 - 

P <x F 



2|C °l ( ^ p)l < r log log A' -f 2(1 2/tt ) log log /,- + 0( I ). 



p=l mod 4 



P 



if log k < b^/YogX. By Merten's theorem, 

^- = loglogX + 0(i; 



p<X 



and thus 



J2^Z < (2/tt — I ) log log .r 4 2(1 - 2/ ~) log log/,- 4 O(l) 



p<X 



p 



□ 
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4. Proof of Theorem [T] 

Here we carry out the steps of the proof as indicated in the introduction. 
First recall that the collision operator can be written 

(16) Q(f,f)(v) =4 / ( T g v (w,9)de) dw, 



where, if we identify u £ S 1 with 9 £ [— n, 7r[, 

g v (w, 9) = q(\w\, oos(e)) (f(v')f(vi) - /(«)/(«.)) , 

and 

1/ = v + w + RqW 
vl = v + w — R s w ; 

as before, w = (u* — u)/2, and -Rg denotes a rotation by an angle 9. Writing 
the Boltzmann equation for two-dimensional velocities, of course we have 
already stepped away from the physically realistic case, but disregarding 
this, a common assumption on q is that 

g(H,cos(6>)) = qi(\w\)q 2 (9) , 

where ?i(|w|) ~ \w\ a for some a £ [0, 1], and where 92(e) ~ |e| -7 for 
some 7 £]1,3[. This corresponds to a molecular interaction by hard in- 
verse power law forces. With the stronger assumption that q\ is smooth 
and strictly positive, it is possible to prove that there is a smooth solution 
f(v , t) to the Boltzmann equation (see ifTUlO . and then this also gives some 
regularity to g(w, 9), in spite of the singularity of q 2 . 

However, much work on the Boltzmann equation has been done with the 
hypothesis that q is bounded or continuous with respect to 9. With that 
assumption, the solution f(v,t) keeps exactly the regularity of the initial 
data. 

Because of this, it is relevant to assume whatever regularity of the so- 
lutions that is needed for the computations. With the aim of making the 
calculations easy, Theorem [T] has been written with unnecessarily strong 
hypothesis. 

To simplify notation a little, let 

/'7T 

G v (w) = / g v (w,9)d9, 



in the continuous case, and for the discrete case (then we assume, of course, 
that v £ hi?) 

G v( h = ^ £ 9v(K,9), 
IC'NICI 
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where 9 is the angle between (' and (. As before, r( |C| 2 ) denotes the number 
of integer points on a sphere with radius \(\. 
Let 

(17) Z KR = {z G Z 2 s.t. |*| < R/h} 

for some R > (this is the most natural example, but other choices might 
be more efficient, as we shall see later). We want to prove that 



(18) 







Q(fJ)(v) - (2hY ^ G»{hQ 

Cez hiR 

when h — > 0, and also make as precise a statement as possible about the 
rate of convergence. 

Theorem 8. Suppose that g v (w, 9) in JToT ) satisfies 

(1) g v (w, 9) is a C 1 -function w.r.t. w 

(2) g v (w, 9) is a C 2 -function w.r.t. 9 

(3) \\g v (-,9)(l + \-\ 2 )\\ L i (dw) <C 

(This holds e.g. if the function f and the crossection q are C 2 .) For given 
R > and h > 0, let Zh,R be as in rti71) . Then given e > there are reals 
R > and h > such that 

Q(fJ)(v) - (2h) 2 J2 G v( h O 

Proof. We still consider Q(f, f) as an iterated integral, and write (for v G 

/iZ 2 ) 

Q(/,/)( V )-(2/ i ) 2 £ G>C) 

= / G v {w)dw-{2h) 2 V G„(/iC) 

(19) + (2h) 2 {G v (h() - G h v (h()) ■ 

From the third part of the hypothesis on g (which is implied by a decay 
of f(v) for large velocities), it follows that for all R > 0, 

(20) / G v (w)dw<^. 

J\w\>R ft 

Continuity of G v (w) would be enough to conclude that 



/ G v {w)dw-{2h) 2 °v( h 

J\W\<R f^rj 



(&z htR 
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when h — > 0. The hypothesis on g v (w, 9) implies that actually G v {-) E C l , 
and there is a constant C 2 such that the difference is smaller than 

(21) C 2 R 2 h = Cmax\d w .G v (w)\ R 2 h. 

w,j 

Next we turn to the difference G v (h() — G^(h(), i.e. of 

( 22 ) ^ [j v{hU)dd ~Hm ^ 2 9v{hceh 

IC'HICI 

(recall that in the second term, 9 is the angle between £' and O- We first 
write the periodic function g v (h(, 9) as a Fourier series, 

g v (h(,e) = Y,9v(Cky k \ 

where 



g v ((,k) = —J g v (h(,9)e-^d9. 

The assumptions on g imply the existence of a constant C 3 so that 

C 3 



(23) 

Then (l2"2l) becomes 



l^(C,fc)l < 



1 + k 2 



ike 



|C'I=|C| IC'NICI 

where the first terms cancel out, and only last sum remains. We next split 
that sum into a part with \k\ < M, and a remainder, which can be made 
small by choosing M large, if g is sufficiently smooth with respect to 9. 
Using 



IC'NICI 



|fc|>M 



< 2— . 

- M 



To find the contribution of this term to ( fl9t . we multiply by (2/i) 2 and sum 
over ( G Zj^ij to find a bound of the form 



(24) 



M 
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For the remaining part, using d23t again, we find a bound of the form 



e lke 



< max 

0<|fc|<M 



s(\C\ 2 ,k) 



KiCP 



EC3 1 
l + Fr(lCl 2 , 
o<|fc|<M vm ' c'ez 2 

IC'HCI 

(25) 

Adding the error terms (l20b . (|2TT) . (l24b and <|25]) gives 

IW,/)(tO-Q fc (/\ /*)(*)! 

(26) 



EC*3 
1 + k 2 

0<\k\<M 



C 3 (2h) 2 max V 

V ^ 0<|fc|<A/ ^ 

Cez hiR 



S(\C\ 2 ,k) 
KICI 2 ) 



In the sum on the right hand side, 

s(|CI 2 ,*0 



r(\C\' 



Yl \ s ( n ' k )\i 

n<(R/h) 2 

and this can be estimated by using Proposition |6] with X = (R/h) 2 . To do 
this, we must require that 

(27) R/h > exp (log(M) 2 /&) 
for some positive constant b. Then there is a constant C 5 such that 

^Ej S (,aoi<c4-je X p(-( 1 --) li - 5r j. 

The last term in (l26l) will always be the dominating one, and at this point, 
it does not give much to try to optimise the choices of R, M and h. Hence 
to achieve an error of magnitude e we 

(1) takei?= y/ACJe, 

(2) observe that we must have h < e/(AR 2 C 2 ) = e 2 /(AC 1 C 2 ), 

(3) choose M = 4R 2 C A /e = QAdC^/e 2 . 

With these choises of R and M, the last term can then be bounded by 

(28) 4C3C5 exp - 1 ) log 7T 

which converges to zero when h — > 0, and so there is an h so small that also 
the last term in (l26b is smaller than e/4. We see that in order to achieve an 

error of maginitude e, one must take h very small: h = o ^exp(— 2 (loge) 2 e^ 2 ^ 1 

(note that d27l) is then satisfied). 

□ 
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5. Some examples and remarks 



From a numerical point of view, the discretization discussed above would 
be far too costly: a discrete velocity model with N velocities would at least 
correspond to a computational cost of 0(N) per time step, because one 
needs to compute a value for each velocity. When the collision term is 
computed by the sum (fT8t . the cost is 0(N 2 ) times some logarithmic factor 
of N (which comes from the summation over the points on the circles). 
And the calculation above showed that N grows exponentially in terms of 
the accuracy, N ~ \ >> exp(e~ c ) for some positive constant c. 

However, rather than estimating the computational cost in terms of the 
number of discretization points used, it is more relevant to give the cost in 
terms of the desired accuracy, given that the discretization points are used 
in an optimal way. The discussion around (fl"8l) suggests that one can reduce 
the computational cost considerably without compromising the order of ac- 
curacy. The poor rate of convergence is due to the approximation of G v (w). 
Generalizing the formula ( fTSt slightly, we can write 



where ph is the local density of Zh- For Z), = {( G Z 2 s.t. \h(\ < R}, 
one has p h = hr 2 . Of course, even more generally one could take a local 
density which is not constant. 

The procedure for constructing a DVM would then be 

• Choose a density p h so that the sum (l2~9l is approximated to the 
desired order. 

• Choose h so small that there exist a set Z h with this density so that 
for all C, E Zh the angular integral is well approximated by the sum. 

For such a model, the computational cost for each velocity would be 
of the order e~ 3 (this estimate is based on the assumption that the cost of 
evalutating the angular integral is e^ 1 , and that the number of velocities 
is 0(e~ 2 )\ lower cost can be acheived if higher order formulas are used 
for approximating the integrals). The problem remains, that a very large 
number of velocities are needed, and hence the total computational cost is 
still excessive. A more challenging task would be to dilute not only the set 
Z h , but to choose in a systematic way subsets Uh C hi? for the discrete 
velocity model, so that Q(f, f)(v) would be well approximated for all v E 
Uh, and to do this in a way that does not require too large tables for storing 
all possible collisions. 

In the last part of this paper, we wish to illustrate the distribution of good 
radii. We then consider ( = (Ci, C2) e Z 2 s.t. < Q (i = 1,2) |C| < 
20000}. This is an extremely large set of points, which corresponds to a 



(29) 
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huge number of velocities (the 0(N 2 ) factor would in this case be of the 
order 10 17 , which is of course absurd) 

Among the circles with radii | £ | in this set, the largest number of points on 
one circle, is 384. In Fig.[H we show all points ( = ( 2 ) with < Q < 
2000, such that the circle passing through ( has more than 72 points. There 
are 36163 points in this set. This is a small fraction of the total number 
of integer points, but they are seemingly well distributed, except near the 
origin. 

Figure El shows points in the range 10000 < Q < 12000. Here the small 
dots denote points on circles having at least 72 points, and the larger dots 
denote points on circles with at least 192 points (there are 141562 and 1 120 
points respectively in these sets). 



2000 
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1600 



1400 
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400 



200 






" -"V- ; ^'"^ ^ 



' 200 400 600 800 1000 1200 1400 1600 1800 2000 

Figure 1 . Lattice points such that circles containing these 
points, contain at least 72 lattice points 
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Figure 2. Lattice points such that circles containing these 
points, contain at least 72 lattice points (small dots), or at 
least 192 points (the larger dots) 
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